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A METRIC FOR GENUS-ZERO SURFACES 


JOEL HASS AND PATRICE KOEHL 


Abstract. We present a new method to compare the shapes of 
genus-zero surfaces. We introduce a measure of mutual stretch¬ 
ing, the symmetric distortion energy, and establish the existence 
of a conformal diffeomorphism between any two genus-zero sur¬ 
faces that minimizes this energy. We then prove that the energies 
of the minimizing diffeomorphisms give a metric on the space of 
genus-zero Riemannian surfaces. This metric and the correspond¬ 
ing optimal diffeomorphisms are shown to have properties that are 
highly desirable for applications. 


1. Introduction 

The problem of comparing the shapes of surface arises in many fields, 
including facial recognition, image processing, brain cortex analysis, 
protein structure analysis and computer vision. It is referred to by 
names such as surface registration, surface warping, best fit, shape anal¬ 
ysis and geometric morphometries. In this paper we introduce a new 
method to compare the shapes of two genus-zero surfaces. The method 
is based on a sequence of two energy minimizations, first minimizing 
the Dirichlet energy to produce a conformal map and then minimiz¬ 
ing a symmetric distortion energy, defined in Section It produces a 
metric dgd on the space of piecewise-smooth surfaces genus-zero Rie¬ 
mannian surfaces, which we call the symmetric distortion metric. In 
addition to giving a distance between any pair of genus-zero surfaces, 
the method also produces an optimal correspondence between them, 
a diffeomorphism whose symmetric distortion energy deffnes the sym¬ 
metric distortion distance. 

A common approach to shape comparison of surfaces in is to 
search for a Euclidean motion moving one surface close to the second, 
and to then measure in some way the setwise difference between the two 
repositioned surfaces. Such approaches are extrinsic, as they consider 
not just the two-dimensional geometry of the surface, but also the 
particular geometric embedding of the shape in space. In extrinsic 
geometry, a hand in different configurations represents very different 
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geometric shapes. From the intrinsic point of view, which we use, a 
hand in different poses represents close to identical geometries. The 
intrinsic approach has significant advantages when comparing surfaces 
that can be flexible. 

Our method of comparing two shapes involves finding an optimal 
diffeomorphism from one to the other, a map that minimizes a mea¬ 
sure of surface distortion. This is often not the case in methods that 
compare surfaces by creating vectors of shape signatures based on fea¬ 
tures such as diameter, curvatures, spectral properties, and spherical 
harmonics. The existence of an explicit correspondence realizing the 
minimal distance is very useful in many applications, and gives an ad¬ 
vantage over methods, both extrinsic and intrinsic, that give distance 
measures without producing surface correspondences. 

A key property of the measure of shape difference that we introduce 
is that it gives a mathematical metric on the space of shapes of genus- 
zero surfaces. A metric on a set A is a distance function d : A x A —?• E 
that satisfies three properties: 

(1) d{S, T) > 0, and d{S, T) = 0 if and only li S = T. 

(2) d(S,T) = d(T,S) 

(3) d(S, W) < d{S, T) + d{T, W) 

These properties are highly desirable for a shape comparison function. 
They imply that the distance between shapes is stable and not overly 
sensitive to noise and measurement error. If S' is close to S and T' 
is close to T, then condition (3) implies that d{S,T) ~ d{S',T'). We 
will introduce a mathematical metric on the space S of genus-zero 
surfaces with piecewise-smooth Riemannian metrics, with two surfaces 
considered equivalent if they are isometric. 

Our method has many additional useful features. It gives a conformal 
diffeomorphism from one surface to the other, useful for applications 
such as texturing. It is well suited to representation of smooth surfaces 
by triangular meshes. The computed correspondence is robust under a 
change of mesh, either from a perturbation of the location of vertices or 
from using a combinatorially distinct mesh. It is intrinsic, so that the 
correspondence between two surfaces does not depend on how they are 
embedded in M^, but only on their Riemannian metrics, and is there¬ 
fore well suited for comparing flexible surfaces that arise in the study of 
non-rigid objects. The method applies to immersed surfaces (surfaces 
with self-intersections) and to surfaces in arbitrary manifolds. Note 
that arbitrary Riemannian surfaces may not be realizable as subsets 
of E^. Furthermore, surfaces whose meshes have intersecting or over¬ 
lapping triangles fit just as well into the framework of our algorithm 
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as embedded surfaces. Pairs of intersecting triangles are common in 
meshes constrncted from point clouds, and are problematic for some 
shape comparison approaches. Finally, we note that onr method can 
be Implemented to be completely automated, and does not rely on 
any preliminary labeling of landmark or feature points. This allows for 
avoidance of errors and costs due to variability of human input. 

Applications of our shape comparison method include: 


(1) Shape retrieval, or hnding nearest fits in an atlas, or database 
of shapes, 

(2) Geometric clustering, 

(3) Alignment of surfaces with different conformations bnt similar 
surface geometry. For example, comparing scans of faces that 
exhibit different facial expressions, 

(4) Alignment of images of one object taken at different times, 
to measnre change over time, and to locate snbregions where 
changes have occurred, 

(5) Statistically sampling surfaces, and averaging to find typical 
surfaces or random surfaces, 

(6) Determining the suitability of a conformal parametrization of 
a surface. Computation of a very large dilation can indicate 
problems in conformal parametrization, resulting in a mesh that 
does not closely align to a modeled underlying surface. 

(7) Transferring a single common mesh to a collection of genns- 
zero surfaces. This in tnrn can be used to interpolate between 
collections of surfaces which are initially described with distinct 
meshes, giving an average shape for a collection of differently 
meshed surfaces, 

(8) Creating a conformal map to use as an initial value or seed 
in shape correspondence methods that allow for non-conformal 
correspondences bnt depend on a good initial correspondence, 

(9) Detecting symmetry. If a surface has reflectional symmetry 
then it and a reflected copy have small distance. Similarly a 
diffeomorphism whose source and image have small distance 
and that is not close to the identity indicates existence of a 
symmetry, 

(10) Coarsening a mesh while retaining surface geometry. The com¬ 
putation of dsci is minimally affected by subdivision or coarsen¬ 
ing of a mesh, so coarsening a mesh will preserve as long as 
the coarsened mesh is geometrically close to the original mesh. 
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We consider here the case where each of the two compared shapes is a 
surface of genus zero, or a topological sphere. The restriction to genus 
zero is appropriate for a wide variety of natural surface comparison 
problems, including facial recognition, alignment and comparison of 
brain cortices, comparing protein surfaces, and geometric identification 
and comparison of objects such as bones and teeth. Note that when a 
comparison is sought between two disk type surfaces, each with a single 
boundary curve, this problem can be transformed into a comparison of 
two spheres. The transformation can be accomplished, for example, by 
gluing a flat disk with appropriate boundary length onto the boundary 
of each of the pair of initial surfaces, turning them into genus-zero 
surfaces. The same idea allows comparison of annuli or more general 
disks with holes. Extensions to surfaces of higher genus can be carried 
out by considering conformal classes of fiat and hyperbolic geometries, 
or by searching for canonical surgeries to reduce a surface to genus zero. 
These will be explored elsewhere. 

Our approach is based on successively minimizing two energies de¬ 
fined on maps between surfaces. We first minimize the Dirichlet energy 
among all maps between the surfaces, yielding a map which is har¬ 
monic. For genus zero surfaces, the harmonic maps exactly coincide 
with the conformal maps, and this step reduces the maps to be consid¬ 
ered from the vast space of all diffeomorphisms to the much smaller, 
but still large, six-dimensional space of conformal maps. We then min¬ 
imize again, this time using a symmetrized energy function introduced 
in this paper that we call the symmetric distortion energy. Minimizing 
this energy amounts to picking an appropriate Mobius transformation, 
as indicated in Figure [T} The symmetric distortion energy gives a mea¬ 
sure of the distance of a conformal map from an isometry. This energy 
is both conceptually natural and efficiently computable, and achieves 
good results in experimental tests. We show in Section]^ that the sym¬ 
metric distortion energy and the metric on surfaces it induces behave 
well as measures of shape similarity. In a related paper we apply the 
symmetric distortion distance to study similarities of shape in biolog¬ 
ical objects such as the surfaces of bones [TH]. Results in that paper 
indicate that this distance is remarkably effective in distinguishing and 
grouping biological shapes. 

1.1. Prior work. Conformal maps from surfaces to the plane have 
become an important tool to visualize and to flatten surfaces, in par¬ 
ticular for surfaces that are topological disks, but also for spheres and 
higher genus surfaces. By mapping a surface region to the plane while 
preserving angles, these methods allow for consistent visualization of 
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Figure 1 . In this example the symmetric distortion dis¬ 
tance between two teeth is calculated. The first surface 
Fi is scanned from the molar of a flying lemur while the 
second F 2 comes from a tree shrew. A conformal map 
minimizing symmetric distortion energy is computed in 
two main stages. In the first stage, conformal maps to 
the round sphere, ci : Fi —>■ and C 2 '■ F 2 ^ S^, 

are computed. In the second step, a Mobius trans¬ 
formation m : —?• is computed that minimizes 

the symmetric distortion energy of the composite map 
/ = C 2 ~^ o m o Cl : Fi ^ F 2 . The energy of / gives the 
symmetric distortion distance between the surfaces. 


locations on highly folded surfaces, and for graphical techniques such 
as texturing. 

Pinkall and Polthier described a mid-edge method of computing dis¬ 
crete conformal maps and applied it to the computation of discrete min¬ 
imal surfaces ng, Bobenko, Pinkall and Springborn gave an approach 
to computing discrete conformal maps based on an energy minimiza¬ 
tion technique |2] . Thurston suggested that discrete approximations of 
conformal maps could be obtained by circle packings. This was carried 
out in work of Hurdal |H] and Stephenson [Til]. Rodin and Sullivan, 
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and He and Schramm established that the limits of discrete maps ob¬ 
tained by circle packings converge to smooth conformal maps [H] , [10] . 
Haker et ah developed a method based on finite elements to compute 
discrete conformal maps [7] . Gu and Yau computed discrete conformal 
parametrizations for surface matching [H]. Jin, Wang, Yau and Gu used 
a stretching energy to create optimal parametrizations of surfaces HU. 
The Gromov-Hausdorff distance has also been used to develop shape 
comparison methods, as in [HI El HE]- 

Recently a series of shape comparison methods introduced that are 
based on optimal transport. Lipman and Funkhouser developed a 
method to find an optimal conformal surface correspondence based 
on a voting scheme that weighs transportation costs ca Boyer et 
ah introduced several innovative methods to compare shapes based on 
minimizing a cost based on optimal transport [Ij. They tested their 
methods on a collection of scanned biological objects, consisting of 
teeth, radius bones, and metatarsal bones from a variety of primates. 
This data was made available and we have used it to test our method 
and to compare its effectiveness to the methods described in [1], and 
to the expert observer data that they provided. 

Earlier work of the two authors used related energies associated to 
conformal diffeomorphisms to compare shapes of brain cortices and of 
protein surfaces na. This earlier work has been improved and further 
developed here. In particular, the optimal diffeomorphism produced 
by the symmetric distortion energy is proven to have values that give 
a metric on the space of shapes, a highly desirable feature not present 
in the previously studied energies. In the discrete setting, the ap¬ 
proach given here has been improved to give mesh independent surface 
comparisons. In contrast, the method given in na required combina- 
torially identical meshes before it could provide a consistent measure 
of distances between shapes, a requirement that restricted the scope 
of applications. These limitations have been overcome in the current 
work. 


2. Aligning smooth surfages 

In this section we develop our method in the context of smooth 
surfaces and mappings. This gives the underlying theory for the sub¬ 
sequent implementations of computational algorithms on triangulated 
or meshed surfaces. 

A Riemannian surface is a smooth 2-dimensional manifold equipped 
with a Riemannian metric, a smoothly varying inner product on the 
tangent space of the surface. An isometry between two Riemannian 


A METRIC FOR GENUS-ZERO SURFACES 


7 


surfaces is a map under which the Riemannian metrics correspond. 
In particular, an isometry preserves angles and distances along the 
surface. Not all angle preserving maps are isometries. Similarities of 
the plane, which stretch all distances uniformly, give an example of a 
non-isometric angle preserving map. The maps that preserve angles at 
each point are called conformal 

Some metric distortion is necessarily present in any construction of 
an alignment from a surface to another when no isometry exists. A 
measure of this distortion is given by the total stretching energy of the 
first snrface as it is deformed over the second. This stretching can be 
measnred by the Dirichlet energy E£){f) of a map / : Fi —)■ F 2 , defined 
by the integral 

fiD(/) = i/ \\df\?dA 

J F-^ 

The maps that minimize Dirichlet energy between two surfaces are 
called harmonic maps. In two dimensions, the class of harmonic dif- 
feomorphisms of spheres, the class of conformal diffeomorphisms, and 
the class of holomorphic (complex differentiable) diffeomorphisms, all 
coincide. We focus here on the angle-preserving property of conformal 
maps. 

A deep result, the Uniformization Theorem, states that a confor¬ 
mal diffeomorphism always exists between two smooth genus-zero Rie¬ 
mannian surfaces Fi and F 2 [I]. However such conformal maps are 
not unique. Each conformal diffeomorphism / : Fi ^ F 2 is part of a 
6-dimensional family. To understand this family we consider the case 
where Fi is the ronnd, radius-one 2-sphere C M^. The space of 
conformal diffeomorphisms from to itself forms the six-dimensional 
gronp PSL{2,C), called the Mobius or linear-fractional transforma¬ 
tions. Any conformal map / : —>■ T 2 can be precomposed with 

a conformal Mobius transformation f ^ to give a new con¬ 
formal map f o (j) : ^ F 2 , and this construction gives the entire 

six-dimensional family of conformal maps from to F2 

A conformal map / : Fi —?• F 2 stretches the metric of Fi at a point 
p E Fi uniformly in all directions. A conformal diffeomorphism then 
defines a real valned fnnction A/ : Fi ^ M+ that measures the stretch¬ 
ing of vectors at each point. The function A/ is called the dilation, and 
is defined by 

f*{92) = A/V, 

where gi, g 2 are the Riemannian metrics on Fi, F 2 respectively, and 
f*{g 2 ) is the metric on Fi obtained by pulling back the metric on F 2 . 
This formula means that for x E Fi, vectors Vi,V2 E FFi in the tangent 
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space of Fi at x, and f*{yi) the image of Vi in the tangent space of F2 
at f{x) under the derivative of /, we have 

fi'2(/=.(Pl),/*(P2))/(x) = ^f‘^{x)gi{Vi,V2)x- 

The conformal factor A/^(x) measures the pointwise expansion or con¬ 
traction of area at x. 

A common measure of the global distortion of a map / : —)• T2 is 
given by the Dirichlet energy of /. This energy is conformally invariant, 
meaning that is preserved by pre-composition with a conformal map. 
The minimal value of the Dirichlet energy over all diffeomorphisms 
between a pair of smooth genus-zero surfaces is obtained when / is a 
conformal diffeomorphism, and in that case the energy is equal to the 
area of F2. For a conformal map with Jacobian determinant |Jac(/)| 
we have ||(i/|p = 2 A/^ = 2 ||Jac(/)|, 

ED{f) = ly [ < df,df > dA > [ |Jac(/)| dA = Area(F2). 

^ Jfi Jfi 

If we take ei, 62 to be an orthonormal frame in a neighborhood of a point 
then the formula for the integrand on that neighborhood simplifies to 

i(ll/.(ei)lt + ll/.fe)in > ll/.(e,)||||/.fe)|| > |Jac(/)| 

Equality holds precisely when / is a conformal diffeomorphism. Since 
the Dirichlet energies of any two conformal maps are equal, we in¬ 
troduce an additional, secondary energy to measure the distance of a 
conformal map from an isometry. This energy emerges naturally from 
the following considerations. If a conformal map / has constant dila¬ 
tion A/ = 1, then df preserves both lengths and angles at each point, 
and thus / is an isometry. So |A/ — 1| indicates the pointwise deviation 
of a conformal map from an isometry. This leads us to the following 
integral that globally measures this deviation, 

f {Xf-lfdA. 

JFi 

We call this the elastic energy of the map /. 

The elastic energy of a conformal map is closely related to the average 
stretching of A/, given by 

EiU) = [ Xf dA. 

Jfi 

Lemma 2.1. A conformal map f minimizes E^ff) among all confor¬ 
mal maps f : Fi ^ F2 if and only if f maximizes Ei among all such 
maps. 
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Proof. For a conformal map / : Fi —>■ F 2 , we have ||4f|P =< df,df >= 
2A/^ = 2|Jac(/)| and ED^f) = Area(F 2 ). Then 


Eelif) 



{Xf - If dAi 



2Xf + 1) dAi 



dAi T 



ldAi-2 



Xf dAi 


ED{f) + Aveei{Ff-2Eff) 
Area(F 2 ) + Area(Fi) — 2Ei(f). 


The first two terms do not depend on the choice of conformal map /. 
Thus Eel is minimized when Ei is maximized. □ 


A somewhat similar functional was introduced by Jin, Wang, Yau 
and Gu m- The integrand in their work has the form {Xf — if rather 
than the {Xf — if used in our definition of elastic energy. They called 
this a “uniformity energy” on conformal maps and applied it to find 
optimal parametrizations of surfaces. However as a result of the fourth 
order term, Lemma |2.1| does not apply to their energy. 

The energy E^i is bounded above in the space of conformal maps 
from El to E 2 by ^(^ 2 ) + A{Ei). The definition of elastic energy can 
be extended to non-conformal maps by taking 


Eelif) = Eoif) + Area(Fi) - [ tr(d/) dAi. 

JFi 

We restrict attention to conformal maps in this paper. 

A drawback of the elastic energy as a measure of shape distortion 
is its lack of symmetry. The elastic energy of / and f~^ may not be 
equal. Nor does the optimality of / imply the same for f~^. We now 
define a symmetrized energy that corrects this shortcoming. 


Definition. The Symmetric Distortion Energy of a conformal diffeo- 
morphism f : Ei ^ E 2 with dilation function A/ : Fi ^ M is 


Esdif) = VEelif) + VEel{f-f 



Note that Egdif) = Egdif ^)- A conformal map that minimizes 
Esdif) gives the optimal correspondence between Ei and E2 that we 
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seek, and the magnitude of Egd for such an optimizing map defines the 
distance between the two surfaces. 

Our next goal is to show the existence of a diffeomorphism that 
minimizes Egdif)- A difficulty is that the pointwise limit of a sequence 
of conformal diffeomorphisms may be discontinuous, or may map all of 
Fi to a single point in F 2 . We need to show that a sequence approaching 
an infimum of the energy does not have this undesired behavior. We 
first examine the special case when / is a conformal map from the 
round sphere to itself. 


Lemma 2.2. If a sequence of Mobius transformations nii : ^ 

has no convergent subsequence, then there is a subsequence for which 
Ei{mi) = 0 . 


Proof. A Mobius transformation is completely determined by the image 
of three points. Let P,Q,R be any there points in S"^. If a sequence 
of Mobius transformations {ruj} takes P,Q, R to points Pi,Qi,Ri, and 
these converge to three distinct points P' ,Q' ,R!, then the sequence of 
Mobius transformations {ruj} is equicontinuous and converges to the 
unique Mobius transformation that takes P,Q,R to P',Q',R'. Thus if 
{nii} has no convergent subsequence then two or more of P,Q,R are 
converging to a single point as z ^ 00 . 

A nontrivial Mobius transformation fixes either one (called para¬ 
bolic) or two (called elliptic or hyperbolic) points on S"^. We divide the 
proof into three cases, according to the limiting behavior of these fixed 
points as i —>■ cxo. 

Case (1): There is a subsequence in which each Mobius transformation 
has a single fixed point. 

Case (2): There is a subsequence in which each Mobius transformation 
has two distinct fixed points, and these converge to two distinct points 
as z ^ 00 . 

Case (3): There is a subsequence in which each Mobius transformation 
has two distinct fixed points, and these converge to a single point as 
i ^ 00 . 

One of these three cases must hold; if Case (1) does not apply, then 
all but finitely many Mobius transformations have two distinct fixed 
points, and by compactness either Case (2) or Case (3) hold for some 
subsequence. In each case we show that the values of Ei on the subse¬ 
quence limit to 0. 

The Riemannian metric of the round sphere, with the north pole 
removed, is isometric under stereographic projection to the plane 
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with the Riemannian metric 


ds^ = 


(1 -|- 

The area form for this metric is 

4 , . 4 


{dx^ + dy^). 


dAi = 


(1 + r^)- 


:dA = 


(1 + r^)' 


-{rdr A dO). 


Case (1): Take a subsequence of consisting of Mobius transforma¬ 
tions with a single fixed point and rotate so that each ruj fixes 
the north pole. Note that conjugating by a rotation does not change 
Ei{mi) and that we still have no converging subsequence since is 
compact. Then each fixes oo in the coordinates given by stere¬ 
ographic projection and thus acts as a translation of the plane. By 
further conjugating by a rotation fixing both 0 and oo, we can arrange 
that that mi{z) = z + Bi, with Rj > 0 € M. Divergence of {m*} implies 
that Bi ^ oo. A computation gives that 

\ 1 1 — ^ 


We first consider the integral of Ei(mi) over the disk Dq = D(0, Bi/2) 
of radius Bj/2 centered at the origin. On Dq we have \z -\- Bi\ > Bi/2, 
so that 


dA\ 


'Do 


'DiO 


1 + 

z\ 

2 

4 


1 + 

^ + Bi\ 

2(1 + 


2^2 


dA 


1 


< 4 


Do 1 + 1^ + 1 + 1^1 

1 


1 


'Do 1 + 1 + 1 ^ 


dA 


dA 


47r 


''Bi/2 


2r 


l + \Bi/2\^J, (l + r2) 

A-tt 

log(l -F Bi^/A). 


dr 


1 + I-Si/2 p 

Note that this integral approaches 0 as S, ^ oo. 

Next we compute over the disk Di = D((—Bi,0), Bi/2) of 

radius Bi/2 centered at (—Sj, 0). The change of variables 2 ; = —w — Bi 
takes Dq to Di and we carry out the corresponding change of variables. 


1 




1 + R + Sdn + R| 


dA = 


'Do 


1 + ItcP l + \w + BA 


dA 


1 
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As before, for w E Dq we have \w + Bi\ > Bi/2, and again this integral 
approaches 0 as S* ^ oo. 

Finally we consider the integral of Ei{mi) over F >2 = C — {Dq U Di}. 
Claim 2.3. \z + -Bj|/| 2 ;| > 1/3 for z G D 2 . 


Proof. If \z\ < 3Bi/2 then 


2 ^ Di 


z Bi ^ Dq 


z + Bi\ ^ I-Si/21 
k| “ |3Si/2| 


1 

3’ 


If \z\ > |3Sj/2| then Bi < 2 |; 2|/3 and \z + Bf > \z\/?> so again \z + 
Bi\/\z\ > 1/3. □ 


Then 


ID 2 


1 T 1^ T SjP 1 + 1^1 


dA < A 


< 4 


'd2 1 + 1 + kl 


dA 


Ic-Do A 

POO 


9 

z^dA 


.00 ^ 

■JBi /2 


= 7271 / — rdr 


14471 


Bi 


2 • 


This term also approaches 0 as S, —?• cx). 

Since Sq U Si U S 2 = C and the integral over each approaches 0 as 
Bi — 7 - 00 , we have shown in Case ( 1 ) that limi_yoo Ei{mi) = 0. 

Case ( 2 ): By conjugating with rotations as before, we can assume that 
each rrii has 00 as an attracting fixed point and also fixes a second 
point Pi- In Case ( 2 ) we assume that hmj^.ooPi 7 ^ 00 , and therefore 
after passing to a subsequence we can assume hmj^.oo Pi = P for some 
point p. We can conjugate rrii by the Mobius transformation 2 ; -E z—pi 
to get a Mobius transformation that fixes 00 and the origin. This gives 
a new sequence of Mobius transformations m' that fixes both the origin 
and 00 and such that aEiinii) < Ei{m!f) < jdEiimi) for fixed positive 
constants a, (3. Thus 


lim Eiijn'i) = 0 lim Eiirrii) = 0 

i^oo i^oo 

and we can assume that rrii fixes the origin and 00 . 

A conformal transformation of the sphere fixing the origin and 00 
with 00 as an attracting fixed point has the form mi{z) = AiZ, |Aj| > 1 . 
Since {rrii} has no convergent subsequence we must have that \Ai\ ^ 
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oo. The value of at is given by 


Xmi 


1 + \z\‘^ 

1 + lApkp 


14 


Ei(mi) 



oo ^ 

(l + ^2)2 

(1 + 4)141 


rdrdO 

4 


1 + 14^4 (1 + 


rdr 


2, r_iiiT_ 

J„ (l + |A|V)(l + r2) 

87r|4|log(|4|) 

I4P-1 ■ 


dr 


Again limj^oo £^i(mi) = 0. 

Case (3): As before, we can first rotate S"^ so that each ruj has the 
north pole as an attracting fixed point. Then mi(z) = AiZ + Bi for 
some Ai,Bi G C with |Aj| > 1. By conjugating by a rotation fixing 0 
and oo we can assume that Bi G M"*". The second fixed point of ruj is 
then Pi — —Bi /(Aj — 1), and since we assume in Case (3) that Pi oo, 
we have limj^oo |4/(4 — 1)| = oo. If Bi is bounded then Aj —)• 1 and 
the sequence rrii = AiZ + Bi has a convergent subsequence, contrary to 
our assumption. Therefore Bi ^ oo and hence also Bj/|Aj| —)■ oo. 

We first compute Ei over the disk Di = D{—Bi/Ai,Bi/\2Ai\). The 
change of variables 2 ; = —w/Ai — Bi/Ai, with Jacobian l/|Ajp, trans¬ 
forms this to an integral over Dq = D{0, Bi/2). Note that for w E Dq 
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we have \w -\- Bi\ > Bi/2. 


fDi 


^rrii (^) dAi 


lAI 


/di 1 + 1 + 

/” |A| 4 


dA 


dA 


< 


CqI + AP 1 + A/A + A/AP IAP 

f I A| _4_ 

Cgl + |t(;P lAP + k + AP 

f I A| _4_ 

'Do 1 + AP IAp + IA/2|" 


dA 


dA 


4| A| 


f27r 



1 


\A\^+\B,/2\^ Jo Jo l + A 

47r|A| : 

|A|2+|A/2|2io l + r^ 


rdrdO 


:dr 


< 


47r|Ai|log(l +A/4) 
|AP + |A/2|2 
167r|Ai|log(l +AV4) 
lAP 


Since |Aj|/A —>■ 0, this quantity approaches 0 as i —>■ oo. 

Over the disk D 2 = D{0, i?j/|2Aj|) we have \AiZ + Bj > Bi/2. 

[ X(mi) dAi = [ - rir~^ —;7T7 — “1 

j D 2 Jd2 1 + IA^ + Ap 1 + l^p 


dA 


< 


lAI 


27r|Aj| rBi/\‘2Ai\ ^ 

T+WWI 

47r|Aj| 


dA 


rdr 


< 


1 + IA/2P 
167r|Ai| log(l + Bf) 


1 + 

log(l + (A/|2A|0 


Bi 


Bi 


This also approaches 0 as i —?• 00 . 

Finally we consider the region = C — {Di U D 2 }. 

Clciiiii 2,4:, In Do wc have jAiZ T BJ ^ Bil2 and jAiZ T Bi I > IV3|- 

Proof. The first inequality follows from z ^ Di AiZ + Bi ^ 

D 2 \AiZ + Bi\ > Bi/2. 
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We consider two cases for the second inequality. If | 2 ;| < ‘iBi/\2Ai 
then z ^ D\ AiZ + Bi ^ Dq and 


\AiZ + Bi 
\z\ 

If \z\ > 3Bi/\2Ai\ then 

2|A^| 


> 


Bi/2 lA 


Bi < 

- 3 

giving the second inequality. 
Then 


3A/|2A 


\AiZ + Bi\ > 


> 


\Aiz\ 


>!£l. 


'Da 


dAi 


lAI 


I Da 


1 + \AiZ + AP 1 + 1^ 


dA 


< 


< 


lA 


'Da IA^ + AP I^P 


dA 


lA 


JDa 

< 24 


|2;/3|(A/2) I^P 

[ ^ 

C-U 2 Bi\z 


dA 


dA 


[ 4 

lBi/\2Ai\ r 


□ 


487r| A| 

Bi 

96|AP 

Bi^ 

Since |A|/A —t 0, this also approaches 0 as i —>■ 00 . 

We have shown in all cases that liruj^oo= 0, proving the 
Lemma. □ 


We now show that Egd achieves a minimum for an appropriate choice 
of conformal map fo '■ Fi ^ F 2 . Let 

I = mi{Egd{f) : f : Fi ^ F 2 is a conformal diffeomorphism} 

A minimizing sequence of conformal maps {/j : Fi —>■ F 2 } is a sequence 
with non-increasing values for Egd such that 


lim Esdifi) = T. 

Z—)-oo 

Theorem 2.5. There exists a conformal diffeomorphism / : Fi —>■ F 2 
with Egdifi) = T. 
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Proof. By the Uniformization Theorem we know there exist conformal 
diffeomorphisms Ci \ Fi ^ S‘^ and C 2 \ F 2 ^ S‘^. The set of all 
conformal diffeomorphisms from Fi to F 2 is given by maps 

f = C2^ omoci : Fi ^ F 2 , 

where m : ^ is a Mobius transformation. Thus we need to 

show that an appropriate choice of Mobius transformation gives a 
minimizer for Egd- 

Let {fi : Fi ^ F 2 } he a minimizing sequence. Then for each /* we 
have Esdifi) < Egdifi) We need to show that {/*} has a subsequence 
converging to a map /o with Egdifo) = F For each i we can write 
fi = cf^omiOCi for some Mobius transformation rrii : ^ For 

the maps fi = cf^rriiCi : Fi F 2 , we have 

E,if,)= f Xf,dA= f 1 dA 

JFi J Fi 

By compactness of Fi, F2 and there are positive constants a, A and 
b,B such that a < < A and b < X^-i < B. Letting c = ab and 

C = AB it follows that 


cEi(fi) < Ei{mi) < CEi{fi). 


In particular, 

Eiifi) —>■ 0 Eiirrii) —>■ 0 Ei{m^ 0 0 . 


Assume now that fi has no conv ergent subsequence 
does rtii or f~^, and Lemma 
Eiiff^) —>■ 0. Recall that 


2.2 


Then neither 
implies that both Ei{fi) 0 and 


Esdif) = VEetif) + VEelif-^) 

= yA(F2) + A(Fi) - 2Fi(/) + yA(F2) + A(Fi) - 2Fi(/-i) 


so that 

Esdifi) 2yArealFy+ArealA). 

Now Esdifi) < 2A/Area(F2) + Area(Fi) — e for some e > 0 so that for 
i sufficiently large, Egdifi) > Egdifi)- This cannot happen for a mini¬ 
mizing sequence, and thus we have a contradiction to the assumption 
that {fi} does not have a convergent subsequence. 

A convergent subsequence of a minimizing sequence gives a new min¬ 
imizing sequence. A convergent sequence of Mobius transformations 
that limit to a Mobius transformation converges smoothly to the limit¬ 
ing map, and as a consequence the conformal maps {fi} also converge 
smoothly to a limiting map /o, whose energy is equal to I. □ 
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2.1. A metric on shape space. We now show that gives a metric 
on the space of spherical shapes S. 

Define a distance function dgd '■ S x S ^ TZ hy taking the minimal 
symmetric energy over the space of conformal maps: 

dsd{Fi, F 2 ) = X = mi{Esd{f)\f : —?• F 2 is a conformal diffeomorphism}. 

Theorem 2.6. The function dgd defines a metric on S. 

Proof. Let Fi,F 2 ,F 3 be three genus-zero smooth surfaces. To show 
that dsd is a metric we need to check that: 

(1) dsd{F^,F 2)>0 

(2) dsd{Fi,F 2 ) = 0 Fi and F 2 are isometric. 

(3) dsd{Fi,F 2 ) = dsd{F 2 , Fi) 

(4) dsdiFi, F 3 ) < dsdiFi, F 2 ) + dad{F2-, Ffij 

The first three properties are direct consequences of Theorem 2.5 and of 
the formula for the symmetric distortion energy Egd- We now establish 
the triangle inequality. Suppose that / : Fi —?• F 2 and g : E 2 ^ E 3 are 
conformal diffeomorphisms with dilations A/ and Xg. Note first that 
g o f : El ^ E 3 is a conformal diffeomorphism with dilation XfXg. We 
establish first the following lemma. 


Lemma 2.7. 

VFelif) + VFelig) > yjEel{go f) 
Proof. Let K = ^/Eel{f) + ^/Eel{g). Then 
= EM) + Eeiig) + 2 ^EM)\/Eei{g) 


— / ~ + / (1 “ 2^3)^ dA 2 -\- 2 a I / (1 - A/)2 dAi / (1 - Ag)- 

JFi J F2 V *^^2 

A change of variables in the second term gives 


/ {l-X^YdA^^ / (1 - A,)Mac(/) dAi = / 

J F2 ^ F\ J F\ 

Therefore, 


K — / (1 ~ Aj) + (1 — A^) Aj dAi + 2 
Jfi 


(1 - XfY dAi Yil - XgYX} dAi. 


The Cauchy-Schwartz inequality can be applied to the second term on 
the right side of the equation above, 


/ (l-\,rdA, / (l-\,fydA,.> / (l-A/)(l-A,)A,<i/l:. 

IFi JFi JFi 


dA 2 . 
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Therefore, 

> / (1 - A/)^ + (1 - dAi + 2 f {1 - A/)(l - Xg)Xf dAi 

JFi JFi 

> / (1 - A/)^ + (1 - Xg)^Xf + 2(1 - A/)(l - Xg)Xf dAi. 

J Fi 

Expansion of this last integrand gives: 

(1 - XfY + (1 - Ag)^Aj + 2(1 - A/)(l - Xg)Xf 
= 1 - 2A/ + Aj + Aj - 2XgX) + X^X) + 2A/ - 2X) - 2XfXg + 2X‘jXg 
= 1 - 2 A/A 3 + A^A^ 

= (1-A/A,)2. 

Therefore, 

> / (1 - XfXgf dAi 

JFi 

> Eeiigof). 

As K and Eei{g o /) are positive, this concludes the proof. □ 


We now return to the proof of Theorem 2.6. From Lemma |2.7 
have: 


we 


VEelif) + VEelig) > VEeligof). 

Similarly, 

VEelif-^) + VEelig-^) > VEelif-^ og-^). 


By summing these two inequalities, we get, 

( 1 ) E^dif) + E,a{g) > E,a{g o f). 

Now let us suppose that fi is a diffeomorphism that realizes dsd{Ei, E 2 ) 
and gi a diffeomorphism that realizes dsd{E 2 , E 3 ), so that dgdiEi, E 2 ) = 
Esdi^fi) <Mid dgd{^E 2 iE^ Egdi^gx)• From Equation 0 we conclude that 

Egdigi o /i) < dsd{Ei, E2) + dsd{E2, E3). 

Since F 3 ) < Esd{giofi), the triangle inequality holds for dgd- □ 


2.2. A scale invariant distance. In many applications it is desirable 
to compare two shapes that are defined only up to scale. Sometimes 
data is presented in units which are unknown and cannot be compared 
to standardized units. It can also be useful to separate scale from other 
aspects of shape comparison. For this reason we introduce a scale in¬ 
variant version of the symmetric distortion metric, which we call the 
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normalized symmetric distortion metric. The computation of the nor¬ 
malized symmetric distance of two surfaces is found by hrst rescal¬ 
ing so that each has surface area equal to one, and then computing 
the symmetric distortion metric as before. This normalized symmet¬ 
ric distortion metric will be used in the mathematical and biological 
measurements that we describe in Section IH 


2.3. Other energies. Arguments similar to those of Lemma [2^ imply 
that many other energy functions on the space of conformal maps also 
realize a minimum value on some explicit conformal map. It is not 
immediately clear whether these energies define a distance metric on 
the space of genus-zero surfaces. 

Definition. The distortion energy of a conformal diffeomorphism 
/ : Fi —7- F 2 with dilation function A = e“ and 1 < p < cx) is given by 


EAf) = 



1/p 


u{x)\^ dAi 


An energy-minimizing conformal diffeomorphism exists for each Ep, 
since if a minimizing sequence has no convergent subsequence then 
A —>■ 0 and |u| —>■ 00 away from small neighborhoods of two points. 
It then follows that Ep{f) —>■ 00 for each p. The contribution of the 
inverse can be added in as before to give a symmetrized energy, which 
is also unbounded outside of a compact set of Mobius transformations. 


Esdpif) = 


r 

i/p 

r 

/ \u{x)\^ dAi 

+ 

/ \u{y)\P dA2 

Jfi 


JF2 


Note that f~^ has dilation 1/A = so that the formula for Ep{f~^) 
appears similar to that of Ep{f), but the function | — m| is evaluated 
on F 2 for the case of f~^. 


3. Meshed Surfaces 

In applications we generally work with surfaces described by meshes, 
or piecewise-flat triangulations, rather then smooth surfaces. These 
metrized triangulations either have coordinates in given for each 
vertex, or have a length given for each edge. In either case a metric 
is determined in which each triangle is flat and each edge has an as¬ 
signed length. The metric is smooth except at the vertices, where the 
surrounding angle may be less than 27r. A conformal map is approx¬ 
imated in this setting by an appropriately defined discrete conformal 
map, as described by Luo ca 

Given a surface of genus zero Fi with a metrized triangulation we 
compute a discrete conformal map f : Fi ^ from Fi to the unit 
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2-sphere in using the algorithm of Bobenko, Pinkall and Spring- 
born j2]. While we have adopted this procedure for the computations 
presented here, we note that other methods of computing discrete con¬ 
formal maps, such as circle packings or the discrete Ricci flow, can also 
be used. 

We consider a triangular mesh A4 = {V,E,T) in a surface. We 
do not restrict its combinatorial type. The geometry of the surface 
represented by A4 is encoded in its edge lengths. A discrete metrie 
on M is a function I defined on the set of edges E of the mesh, which 
assigns to each edge Cij a length lij so that the triangle inequalities are 
satisfied for all triangles in T. 

When working with meshes and discrete metrics, the elastic energy 
integral is approximated by a sum over the mesh. We consider two 
triangular meshes A4.\ and Ad 2 in with possibly different combina¬ 
torics and different geometries. The geometries are encoded either in 
the positions of the vertices or an assignment of lengths to the edges. 
Given a transformation / : Fi —?• F 2 , in [12] we worked with an elastic 
energy given by 

where indicates the length of the edge in Fi, /(/(ejj)) the 

length of the image of this edge in F 2 , and the sum is over all edges of 
the mesh on Fi. When / is conformal and the mesh is close to uniform 
then the quantity l{f{eij))/l{eij) approximates A/, and this sum is an 
approximation of the symmetric distortion energy. However the sum 
is dependent on the size of the mesh, increasing with the number of 
edges. To make this quantity mesh independent, we weight the terms 
bythe area of the region to which each edge contributes. This leads to 
the following formula for a mesh independent elastic energy. 


Hf) = E 

Cij 





Here Aij is the sum of the areas of the two triangles adjacent to edge eij 
and the sum is over all edges of the mesh on Fi. The weighting factor 
(1/3) assigns to each edge the portion of the area of the two adjacent 
triangles obtained by dividing the triangles into three pieces. 

The symmetric distortion energy in the discrete setting is then ob¬ 
tained by summing over the edges E of the Fi mesh and the edges E' 
of the mesh on F^. 
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3.1. Procedure and Implementation. We begin with two combi¬ 
natorial surfaces Fi,F 2 with metrized triangulations ti,T 2 . We then 
implement the following steps. The process is indicated in Figure 

(1) Construct conformal maps to the unit sphere. 

We use the methods of [2] to construct discrete conformal maps 
Cl : Fi ^ and C 2 '■ F 2 ^ from each of a pair of genus-zero 
surfaces Fi, F 2 to the unit sphere S"^. 

(2) Move the centers of mass of the vertices to the origin. 

This step is done for numerical stability. We compose Ci with a 
Mobius transformation mi and C 2 with a Mobius transformation 
m 2 so that the vertices of mioci(Ti) and m 2 0 C 2 {T 2 ) have centers 
of mass at the origin. This step is done to prevent a choice of a 
conformal map Ci which pushes most of the vertices into a small 
neighborhood of one point on the sphere. Any choice of confor¬ 
mal map from Fi —>■ and from F 2 is theoretically valid 

for our method, but some are computationally problematic. 

(3) Map the source mesh onto the target surface. 

A Mobius transformation m : —?• induces a map C 2 ~^omo 
Cl of the vertices of ti to F 2 . Given a vertex Vj in Fi, we identify 
its image n' in the spherical mesh Ci(Fi). We then locate its 
image n" = on the spherical mesh 02 (^ 2 ) and transfer this 

point to the surface F 2 by applying The image of a point 
that is not a vertex is specified using barycentric coordinates of 
the simplex that contains the point. 

(4) Find an optimal Mobius transformation. 

We search for the Mobius transformation m : —>■ that 

gives rise to a closest to isometric mapping among conformal 
maps between the two surfaces of interest, by searching for a 
global minimum of Egd as given in Equation (|^. We obtain a 
candidate as the solution of a non-linear optimization problem, 
via a steepest descent approach to solve this problem. Steepest 
descent methods are generally fast, but sensitive to local minima 
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and thus dependent on the choice of an initial approximation 
to a solution. 

A random or fixed initial guess, such as the identity transfor¬ 
mation, is likely to lead to a non-optimal local minimum. Each 
initial guess is determined by specifying the images of three 
fixed points on Fi. We can get a collection of initial assignments 
comparable in density to the size n of the mesh on F 2 by choos¬ 
ing all possible assignments for these three points that send 
them to vertices of F 2 . The set of possible choices is then 0{n^), 
which is prohibitive for large meshes. We use a procedure devel¬ 
oped in [12] to automatically generate a collection of reasonable 
initial starting points. The method uses ellipsoid approxima¬ 
tions to Fi and F 2 to give the initial alignment. Each ellipsoid 
approximation generates six points on each surface, correspond¬ 
ing to extremal points where the three coordinate-axis meet the 
surface. We label these points and . An 

initial choice of Mobius transformation is uniquely determined 
by the image of three of the points on F\. We have six choices of 
where to initially map x\, namely any of The point 

x\ is then assigned to the antipodal point on the F 2 ellipsoid. 
We then have four choices of where to map namely to any of 
the four points orthogonal to the image of the first point on the 
ellipsoid. The image of z\. is then determined by orientation. 
Thus we have a total of 24 choices of initial mappings that are 
orientation preserving. Once the image of three points is speci¬ 
fied, a unique Mobius transformation is determined, and this is 
used as one of our collection of initial maps. 

If we also want to consider orientation reversing correspon¬ 
dences, then we first reverse the sign of each z-coordinate of Fi 
and then reapply the process using the reflected surface. This 
gives a total of up to 48 initial correspondences in the unori¬ 
ented case. 

We then apply steepest descent based on the symmetric dis¬ 
tortion energy to find an optimal conformal transformation. We 
use Equation Q to compute the symmetric distortion energy of 
/ and and the gradient of this energy to find a minimum 
value. We compute the symmetric distortion distance as the 
smallest value found for the symmetric distortion energy, and 
use the associated diffeomorphism as an approximation of the 
symmetric distortion energy minimizing map. 
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4. Geometric Computations 

In this section we explore the geometric meaning of the h^d-distance 
by computing it for pairs of well-understood geometric objects. This 
allows us to develop a sense of what dgd is measuring. We measure 
the dsrf-distance between spheres of varyiing radii, ellipsoids of vary¬ 
ing principal axes, and surfaces of varying roughness. We also study 
the effects of decreasing the density of a mesh and of changing the 
orientation of a surface. 

Features that we would like to see hold for dgd to allow its use as a 
robust shape measurement tool are: 

(1) High sensitivity to small changes in area, 

(2) High sensitivity to small changes in shape, 

(3) Low sensitivity to small amounts of noise, 

(4) Mesh independence, 

(5) Low sensitivity to deformations that preserve intrinsic surface 
geometry, 

(6) Ability to distinguish an object from its reflection. 

We show by a series of computational experiments that dgd exhibits 
highly favorable behavior for each of these features. 

Feature (1), sensitivity to area change, can be valuable in some set¬ 
tings, such as measuring the growth of an organism or of a tumor over 
time. In other settings we want to consider only shapes up to scale, such 
as when scans are obtained without a consistent measurement scale. If 
we want to ignore the effect of changing area, we can normalize all areas 
to one by rescaling. Feature (2), sensitivity to small changes in shape, 
can be measured in a variety of ways. Below we investigate the effect 
on dsd-distance of the deformation of a sphere to an ellipsoid which is 
stretched along one axis while maintaining constant area. Feature (3) is 
essential for robust distance measurements that are not unduly affected 
by small amounts to noise or measurement error. Nose sensitivity is 
measured by looking at the effect on d^^-distance of random perturba¬ 
tions of the vertices of a sphere. The mesh independence property of 
Feature (4) implies that ds^-distances are not dependent on the choice 
of a mesh or triangulation used to represent a surface. This allows for 
comparing the geometric similarity of objects having meshes of varying 
density and combinatorial type, subject only to the mesh accurately 
representing the surface. Feature (5) is important for the comparing 
of flexible surfaces, such as the surfaces of proteins, faces and animals 
that take on different configurations or poses. Feature (6) allows for 
the d^rf-distance to distinguish objects that differ only in chirality, such 
as left and right hands, or left and right molars. 
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4.1. Area rescaling. In many applications shapes are presented with¬ 
out scales. For example, two medical images produced with different 
machines can describe the same shape in different coordinates whose 
relative magnitudes is not known. Thus it is often convenient to hrst 
rescale each of the two surfaces being compared so that they have the 
same area, which we can take to be equal to one. 

However in some cases it is useful to measure the effect of a change 
of scale. For example one may want to measure the growth of an object 
over time. When scale is the only difference between two shapes Fi and 
F 2 then Egd measures an integral of the stretching required to enlarge 
one to ht the other. The formula for the energy required to perform 
such ra rescaling can be directly computed. The distance dsd{Si, S 2 ) be¬ 
tween spheres Si of area Ai and S 2 of area A 2 whose optimal alignment 
is realized by rescaling can be computed using Equation Q, giving 


dsd{Sl,S2) — 2 \fA2— 


4.2. Area preserving shape deformation. To measure the effect 
of global changes in shape on the distance dgd between two surfaces of 
equal area, we ran a computation that measured the distance between 
surfaces in a family of ellipsoids from the unit sphere in Two of 
the principal axes of each ellipse are held hxed at radius one, while 
the third is varied from 1/100 to 10. The areas of all surfaces are 
then normalized to one by appropriately rescaling the surfaces, and the 
minimal symmetric distortion energy is then computed. The results are 
indicated in Figure where the distance of each ellipsoid from the unit 
sphere is given as a function of the length of the third axis. Note that 
the distance increases linearly near the point where both surfaces are 
unit spheres, indicating that dgd has the ability to differentiate small 
changes in shape when the two surfaces are close to isometric. This 
feature is highly desirable for the use of dgd as a tool for classifying 
surfaces, as it shows that near-similar surfaces can be differentiated. 
In contrast, the sphericity, a common measure of similarity to a round 
sphere that compares the isoperimetric ratio of a surface to that of a 
sphere, is insensitive to small changes in shape near an isometry, as 
shown in Figure]^ 

A limitation of the current implementation of our computation of 
the dsrf-distance is visible in this experiment. Our method of discrete 
approximation involves measuring the effects of stretching edges of a 
mesh, and this leads to maps that try to avoid sending edges of the 
mesh far out into spikes or protrusions. These issues occur in ellipsoids 
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with one principal axis stretched by a factor close to 10, as shown in 
the graph of Figure]^ (C). 



Figure 2. The effect of global shape on the distance 
between surfaces is indicated by deforming one axis of 
an ellipsoid and computing the distance to a unit sphere. 

All surfaces are scaled to have area one, so it is only 
shape differences that are being measured. For ellipses 
that are close to the round sphere to which they are 
being compared, the dsd distance shown in (C) is more 
sensitive to small changes than the sphericity shown in 
(B). Computational issues arise for the meshes used here 
when A becomes close to 10, as seen in (C). 

4.3. The effect of noise. In reconstructing surfaces from scanned 
data, one often encounters errors in the location of vertices on a surface. 
These variations of vertex positions are local in nature and do not affect 
the overall shape of a surface, but can cause crinkling and spiking effects 
locally. To measure the effect on such noise related local deformations 
of a surface, we added Gaussian noise to the surface of a sphere and 
measured the dgd distance of the resulting surface from a round sphere. 
The mesh used to represent the unit sphere had mean length 0.008 and 
we added Gaussian random radial noise to each vertex, with standard 
deviation equal to a multiple N of the average mesh edge length. The 
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results, shown in Figure are extremely promising. They indicate that 
a random perturbation whose standard deviation is between zero and 
the average edge length of the mesh is recognized by dgd as being close 
to a round sphere. 



Figure 3. The effect of noise, or local deformation of 
the geometry, on the distance between two surfaces is 
indicated by adding Gaussian random noise to the ver¬ 
tices of a unit sphere. Again all surfaces are scaled to 
have area one. N indicates the standard deviation of the 
Gaussian as a multiple of the mean edge length, equal to 
0.008 in this example. 


4.4. Subdividing a mesh. To understand the effect of the choice of 
mesh on dsd-, we experimented with the effects of simplifying a mesh 
by removing points. We take for our first surface Fi a sphere S\ whose 
surface is discretized with 1000 points, distributed uniformly on the 
surface. We compare Fi with a series of spheres having varying numbers 
of mesh points. All surfaces are scaled to have area one 

In Experiment 1 a second sphere is represented with N vertex points 
on its surface that are placed randomly, for values of A” up to A” = 1000. 
Larger values of A, up to A = 16000, are obtained by subdividing 
each of the 1000 triangles into either four or 16 similar triangles. The 
experiment is repeated 50 times for each value of A and the average 
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Figure 4. Distances of spheres represented with vary¬ 
ing numbers of mesh points from the unit sphere with 
1000 uniformly distributed vertices. The upper curve 
gives the dsd distance as a function of the numbers of 
randomly distributed vertices. The lower curve gives the 
dsd distance as a function of the number of uniformly dis¬ 
tributed points. Note that as the number of uniformly 
distributed points drops from 10,000 to 1,000, the dgd- 
distance remains near zero, indicating that a change in 
mesh does not affect the distance when there are enough 
points to accurately model the surface. For low num¬ 
ber of vertices, the graphs reflect a larger deviation from 
roundness of the surfaces represented by the mesh. 

and standard deviation of dgd is obtained for these 50 samples. The 
number of vertices N varies from 10 to 16000, indicated by the upper 
plot in Figure]^ In Experiment 2 the second sphere is again represented 
with N points on its surface with N between 10 and 16000, but this 
time the positions of the points are optimized to give a distribution 
that is as uniform as possible. The resulting distances are shown by 
the lower plot on the Figure]^ The results in both cases indicate that 
a change of mesh does not affect the dg^-distance as long as enough 
vertices are kept to maintain a close approximation of the underlying 
geometrical surface. For uniformly distributed points, the number of 
points required to densely approximate the surface of the sphere is 
smaller than for randomly distributed points, causing a gap between 
the two graphs. The location of the vertices of the meshes has no effect 
once there are enough to accurately capture the geometry of the round 
sphere. 

4.5. Chirality and Reflections. The ds^-distance measures the sym¬ 
metric distortion energy of an orientation preserving diffeomorphism. 
This distance can be reduced significantly if we also allow orientation 
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reversing diffeomorphisms. For example, comparing a right hand and 
mirror-image left hand with dgd will give a non-zero distance. There are 
circumstances when we want to ignore this difference in orientation, or 
chirality. 

We can specify that we wish to incorporate into our shape analy¬ 
sis either only orientation preserving diffeomorphisms, or alternately 
both orientation preserving and orientation reversing diffeomorphisms. 
To allow for orientation reversing correspondences, when comparing a 
surface Fi to F 2 , we add an additional surface Fi which we also com¬ 
pare to F 2 . The surface Fi is obtained by reflecting Fi, computed by 
multiplying the ; 2 -coordinate of each vertex of Fi by —1. This gives 
twice as many candidates for an Egd minimizing map, and may lead 
to a smaller distance. We denote the distance of two surfaces given by 
minimizing in this larger class of potential correspondences by dgd, so 
that dsdiFi, F2) = mm{dsdiFi, Fa), Fa)}. 

To see the effect of adding orientation reversing diffeomorphisms, 
we model a right hand by a surface Fi which is a sphere with three 
protrusions, in the direction of the i,j and k vectors. For Fa we take 
a sequence of surfaces where the protrusion in the direction of j is 
rotated in the xy-plane through — z to —j. Its final position represents 
a surface isometric to the reflection of Fi. We compute the distances 
dsd from Fi to each surface in this family, and then the distances dgd 
from Fi which allow for orientation reversal. The results are shown in 
Figure All surfaces are scaled to have area one. 

An interesting example of this phenomenon occurred in an analysis 
of a collection of teeth taken from a variety of primates, both simians 
and prosimians. The high effectiveness of dgd in measuring similarities 
and differences between such biological shapes is described in m A 
typical set of dgd distances between teeth from the same and from 
different families, is shown in Figure 

Data describing the geometry of a collection of teeth was obtained 
from the study of [1] and we are grateful to Y. Lipman for making it 
available to us. The data contained both left and right teeth. The 
distance between two teeth can be highly affected by the choice of 
whether to allowing orientation reversing correspondences, as indicated 
in Figure 

When only orientable alignments were allowed, the ds^-distance was 
not as effective as either a human observer or as the continuous Pro¬ 
crustes distance described in [Ij at discriminating between the teeth of 
simians and prosimians. The effect of allowing both orientation pre¬ 
serving and reversing maps is seen in the ROC analysis in Figure]^ In 
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Figure 5. Distances of spheres with three protruding 
bumps of heights 1.2 (along the positive y-sods), 1.4 
(along the positive x-axis), and 1.6 (along the positive 
2 ;-axis). The bump of height 1.2, initially facing straight 
out along the y-axis, is rotated through the negative x- 
axis around to the negative y-axis. In (A) the resulting 
dsd-distances are graphed with only orientation preserv¬ 
ing diffeomorphisms allowed. The d^d-distance increases 
initially, but then drops slightly as the bump keeps rotat¬ 
ing to the opposite side of the sphere. In (B) orientation 
reversing diffeomorphisms are also allowed, and the dgci- 
distance drops down to 0 after a rotation of tt. In (C) 
and (D), shaded red areas indicate areas of larger stretch¬ 
ing or compression of the domain, shown on the image 
surface. 


this statistical test, the effectiveness of a distance at predicting mem¬ 
bership in a common family is given by the area under a curve, with 
greater area indicating higher effectiveness. The ds^-distance measured 
only with orientation preserving alignments was not as effective as other 
methods at correctly identifying teeth from the same family (dashed 
red curve in Figure [^. This occurred because both left and right 
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Flying lemur B (0.26) 


Tree shrew B 


(0.55) 


(0.55) 


(0.25) 


Figure 6. A computation of d^^-distances between four 
teeth. Two of the teeth are from flying lemurs and two 
are from tree shrews. All surfaces are scaled to have area 
one. The distances between teeth from the same fam¬ 
ilies are lower than between teeth of different families. 
The shaded red areas indicate areas where there is large 
stretching or compression on the domain (surface A), in¬ 
dicated on the image (surface B). Note that distances 
between pairs of teeth from the same family are smaller. 


AlO 


A13 


Q06 


AlOR 


mmmm 


AlO 0.0 
AlOR 0.36 


0.38 

0.18 


0.30 

0.22 


0.36 

0.0 


Figure 7. All teeth here belong to euprimates, with 
AlO and A13 in one family and Q06 in a second. The 
dg^-distance between molar AlO and molar A13 is 0.38. 
Despite belonging to different families, the distance be¬ 
tween molar AlO and molar Q06 has the smaller value of 
0.30. This seeming mismatch is resolved by considering 
chirality. While AlO and A13 are from the same fam¬ 
ily, they have different orientation, or handedness, while 
AlO and Q06 share the same orientation. The reflected 
tooth AlOR has smaller d^^-distance to A13 (0.18) then 
to Q06 (0.22), indicating that d^^-distance is capturing 
information about the family to which the tooth belongs. 
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molars were included in the data set. When orientation reversing dif- 
feomorphisms were also allowed, the Js^-distance performed as well as 
the other methods (solid red curve in Figure [^. The results indicate 
that geometric differences between left and right molars within the 
same family can be larger than those between right molars from two 
different families. 


Genus 



Figure 8. A statistical test of the effectiveness of a dis¬ 
tance in identifying members of two subgroups is given 
by the area under a curve in a ROC analysis. Greater 
area indicates more effectiveness. The solid red curve 
results from dgd, with orientation reversing diffeomor- 
phisms allowed in comparing shapes, while the dashed 
red curve restricts the computation of dgd, to orientation 
preserving maps. The dashed curve resulted because dgd 
distinguished left and right molars from the same family 
that were highly similar after reflection. 


5. Conclusions 

We have described a new method of comparing the shapes of two 
Riemannian surfaces of genus zero. We introduced the notion of sym¬ 
metric distortion energy and established the existence of a conformal 
diffeomorphism between any pair of genus-zero surfaces that minimi zes 
this energy among all conformal maps. We then established that the 
value of the symmetric distortion energy on the minimizing map leads 
to a metric on the space of shapes. We described how to implement 
this method and there results of experiments performed with such an 
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implementation. These experiments indicate that the symmetric dis¬ 
tortion energy has properties that are highly desirable for many classes 
of applications. 
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